Universal Finite Size Scaling Functions in the 3D Ising Spin Glass 



Matteo Palassini 

Department of Physics, University of California, Santa Cruz, California 95064 
Scuola Normale Superiore and INFM, 56100 Pisa, Italy 

Sergio Caracciolo 
Scuola Normale Superiore and INFN, 56100 Pisa, Italy 
(Submitted February 10, 1999) 

We study the three-dimensional Edwards- Anderson model with binary interactions by Monte Carlo 
simulations. Direct evidence of finite-size scaling is provided, and the universal finite-size scaling 
functions are determined. Monte Carlo data are extrapolated to infinite volume with an iterative 
procedure up to correlation lengths ^ ~ 140. The infinite volume data are consistent with a conven- 
tional power law singularity at finite temperature Tc. Taking into account corrections to scaling, we 
find Tc = 1.156 ± 0.015, = 1.8 ± 0.2 and 77 = -0.26 ± 0.04. The data are also consistent with an 
exponential singularity at finite Tc, but not with an exponential singularity at zero temperature. 

PACS numbers; 75.10.Nr, 64.60.Fr, 75.40.Mg, 75.50.Lk 



The critical properties of the Ising spin glass in three 
dimensions are still not very well understood. Numerical 
simulations have led to some progress [QJ^ , but have been 
hampered by technical difficulties. Large-scale Monte 
Carlo (MC) simulations at correlation length <^ « 10 lat- 
tice units 1^-^ are consistent with both a continuous 
phase transition with power-law divergence of ^ at fi- 
nite temperature T — T^, and an exponential divergence 
at T = 0, which is expected at the lower critical dimen- 
sion. High-statistics MC simulations of smaller systems 
[^"H give a certain evidence of a Tc 7^ transition with 
an ordered spin glass phase below Tc , but cannot exclude 
neither an exponential divergence at T = 0, nor a line of 
critical points at T <Tc ^ (^[H^ , as in the Kosterlitz- 
Thouless theory of the 2D XY model. Understanding 
whether an ordered spin glass phase exists in three di- 
mensions is clearly an issue of major interest. 

In this work, we study the 3D Ising spin glass with 
an approach, based on finite-size scaling (FSS) and MC 
simulations in the paramagnetic phase, introduced in 
Ref. [|lO| (see Ref. for similar methods) and so far 
applied to non-disordered systems. Let us summarize 
our main results, (i) We provide a direct test of the FSS 
hypothesis, independent of the nature of the divergence 
in the infinite system. In particular we determine, for the 
first time to our knowledge, the universal FSS functions, 
(ii) We demonstrate the effectiveness of an iterative pro- 
cedure to extrapolate the MC data to infinite volume, 
that allows us to reach ^ « 140. (iii) Exploiting the 
higher range of ^, we show that an exponential diver- 
gence at T = is excluded, but we still cannot decide 
between a power-law divergence at Tc 7^ and a line 
of critical points terminating at Tc ^ 0. (iv) Under the 
hypothesis of power-law divergence, we show that cor- 
rections to scaling are important and we estimate Tc and 
the critical exponents. 



Model and FSS method — We consider the 3D 
Edwards-Anderson model, whose Hamiltonian is 



(1) 



where ax are Ising spins on a simple cubic lattice of linear 
size L with periodic boundaries, and Jxy are independent 
random interactions taking the values ±1 with probabi- 
lity i. The sum runs over pairs of nearest neighbor sites. 

Let ^(T, L) be a suitably defined finite-volume correla- 
tion length, and let 0{T,L) be any singular observable, 
such as ^(T, L) itself or the spin-glass susceptibility (see 
below). Then FSS theory jTj predicts that 



OiT,L) 
C(T, 00) 
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where fo is a universal function and corrections to FSS 
are neglected. From Eq. (||) one obtains the relation 



0(T, 2L) 
0{T,L) 



- FoU{T,L)/L 
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where Fq is another universal function and only finite- 
volume observables are involved. Our approach works 
as follows (see Ref. jl^ for details). We make MC 
runs at numerous pairs (T, L), (T, 2L) and we plot 
0{T,2L)/0{T,L) versus ^{T,L)/L. If all these points 
fall with good accuracy on a single curve — thus verifying 
the Ansatz (||) — we choose a smooth fitting function Fo . 
Then, using the functions F^ and Fo , we extrapolate the 
pair (^, O) iteratively from L ^ 2L ^ 2^L ^ . . . — ^ 00. 

Computational details — We simulate the model in 
Eq. (|l]) with the heath-bath algorithm. We measure Qx = 
axTx and q = L^^ qx from two independent replicas 
((7, r) with the same Jxy We choose as a definition of 
^(T, L) the second-moment correlation length 
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[S{0)/S{p) - 1] 
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2sin(|p|/2) 
where S{k) is the Fourier transform 

S{k) 



(4) 



(5) 



(arguments T, L are omitted) and p = (0, 0, 2tt/L) is the 
smallest non-zero wave vector The spin-glass sus- 
ceptibility is xsg{T, L) EE L^{q^) = 5(0). The symbol (•) 
represents a double average over thermal noise and J^y^ 
which is estimated from Ng samples with different Jxy 

The runs are done on a Cray T3E parallel computer 
with a fast code that exploits the parallelism of spin glass 
simulations. The binary variables ax and Jxy at corre- 
sponding sites of 64 samples (each represented by a single 
bit) are stored in a 64-bit integer variable, and 64 ax^s 
are updated simultaneously with only 31 logical instruc- 
tions and one random number 0. Average speed on a 
single processor (PE) is 4.5 x 10^ spin updates per second 
(DEC Alpha EV5, 600 MHz). The PEs are arranged in a 
virtual parallelepiped along whose axis we can distribute 
independent groups of 64 samples, different "slices" of a 
large lattice, and different temperatures. We typically 
used 32 to 128 PEs. Equilibration of the runs is verified 
with the criterion introduced in Ref. ]l5[ |. The sizes sim- 
ulated range from L = 4 to L = 48, from which we form 
104 pairs {T,L), {T,2L). In Table I some parameters 
of the simulations are given. The equivalent of about 2 
years of computer time on a single PE was employed. 

TABLE I. Maximum number of samples Ns, minimum 
temperature Tm and Monte Carlo sweeps (MCS) performed 
at Tm , as a function of the size L. 



L 


4 - 8 


10 


12 


16 


24 


32 


48 




1920 


1536 


960 


448 


448 


448 


64 




.9401 


.9793 


1.0936 


1.1642 


1.2059 


1.3397 


1.4084 


MCS/IO"^ 


3 


10 


10 


10 


10 


10 


10 



FSS analysis — In Fig. |^ we show that, within our 
statistical accuracy, the FSS Ansatz (||) is well verified for 

— XSG and O — No systematic deviations from the 
curves are detectable, but data at L = 4, not displayed in 
Fig. 1^, are significantly outside the curves ior ^{T,L)/L = 
X > 0.2. We verified that other observables, such as 
the Binder ratio, also satisfy Eq.(^). We emphasize that 
FSS was not assumed a priori and that Eq.(H) contains 
no adjustable parameters. Furthermore, no particular 
dependence of ^ and XSG on T was assumed. 

We fit the data in Fig. |l| to two functions F^^^ , of 
the form F{x) = l + J2i=i n exp(— i/a;), obtaining good 
fits with n — 3 or 4 (goodness of fit parameter Q > 0.9). 
Using f^sG' ^6 then compute xsciT, oo), ^(T, oo) 
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FIG. 1. Finite-size-scaling plot with the form in Eq.(|^ for 
a.) O = XSG and b) O = ^. Error bars (estimated with a 
jackknife procedure) are one standard deviation. 

with the iterative procedure described above. In Table || 
we show that extrapolations from different L are consis- 
tent, providing a test of the method. In our final analysis, 
we take the weighted average of the extrapolations from 
different L. An implicit assumption of the iterative pro- 
cedure is that the Ansatz (^ with a given function Fq 
will continue to hold as L ^ oo. This assumption could 
fail if the system exhibits a crossover at large L, as in 
any FSS analysis. However, as shown in Table ^ ex- 
trapolations at T = 1.4084 from small L are consistent 
with data from large L, which have little or no finite-size 
effects. We therefore believe that a crossover is unlikely. 

In order to test for systematic errors due to corrections 
to FSS, we repeated the analysis excluding L = 5, 6 from 
the fits of -F^sG ; ^? and we found that extrapolated data 
change within their error bars. We have a good control 
on the extrapolated data up to ^ « 140; at lower tem- 
peratures the statistical errors become quite large, and 
the data are more sensitive to the region of high x, where 
there are few data from large L. (The largest x used for 
the extrapolations is a; = 0.57, from T — 1.2059, L — 5). 

In Fig. we show that with our extrapolated data 
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Eq.(^ is satisfied remarkably well, providing a further 
test of the method. If O ~ S^'^l'' as ^ ^ oo, then 
foi^) in Eq.(||) must satisfy fo{x) ~ x~'^°/'^ as x ^ oo. 
As shown in Fig. |^ (insets), our curves indeed have a 
power-law asymptotic decay, with negative slopes ^jv — 
2.30 ± 0.08 in Fig. |(a) and « 1 in Fig. |(b). 

We emphasize the universality of the scaling functions 
in Fig. |l] and ||. It would be interesting to determine the 
same functions for different distributions of the Jxy, in 
order to test for possible violations of universality ||l6|,|l^ . 



TABLE II. Examples of measured and extrapolated values 
of the correlation length and the spin-glass susceptibility. See 
Ref. llOl] for how to estimate error bars of extrapolated values. 



r 


L 




C(T,oo) 


Xsg{T,L) 


Xsg{T, oo) 


1.2059 


5 


2.85(7) 


120(60) 


36.1(3) 


1.8(5) xlO'' 




6 


3.42(6) 


150(60) 


55.1(5) 


2.8(8) xlO^ 




8 


4.47(6) 


126(30) 


103(1) 


1.9(5) xlO^ 




10 


5.57(6) 


146(30) 


171(2) 


2.8(8) xlO^ 




12 


6.60(8) 


143(30) 


260(3) 


2.7(8) xlO^ 




16 


8.60(15) 


131(30) 


473(10) 


2.1(8) xlO^ 


1.4084 


5 


2.28(5) 


8.5(9) 


25.4(2) 


4.3(3) xlO^ 




6 


2.66(4) 


8.7(7) 


36.1(4) 


4.6(3) xlO^ 
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3.33(4) 


8.4(4) 


60.0(6) 


4.3(2) xlO^ 




10 


3.94(4) 


8.4(3) 


88.3(1) 


4.3(2) xlO^ 




12 


4.51(5) 


8.6(3) 


120(2) 


4.5(2) xlO^ 




16 


5.46(8) 


8.6(3) 


178(4) 


4.4(2) xlO^ 




24 


6.60(11) 


8.1(2) 


269(6) 


4.2(2) xlO^ 




32 


7.16(15) 


7.8(2) 


320(9) 


3.8(2) xlO^ 




48 


8.6(6) 


8.9(7) 


404(30) 


4.3(3) xlO^ 



Nature of the phase transition — We now compare 
our extrapolated data with the following scenarios: (i) 
a Tc / continuous phase transition; (ii) a line of crit- 
ical points terminating at Tc ^ 0, with an exponential 
divergence as T — > r+; (iii) an exponential divergence 
at T = 0. The last two scenarios imply a lower critical 
dimension exactly equal to three, 
(i) We fit our data to 

e(T) = (T - TJ-'^ [1 + ae (T - T,f] (6) 
XsaiO^be-'' [l + dr^] (7) 

with fixed correction-to-scaling exponents and A . 
In the fit we include data with f > ^m, varying in 
order to test the stability of the fits. Without the cor- 
rections to scaling (a^ = d = 0), the quality of fits is 
good for > 3 — 4 (Q w 1), but fit parameters (notice- 
ably Tc, and rj) show small systematic variations with 
in the whole range available. Including the correc- 
tions, we obtain excellent and stable fits with I < < 2 
and 1 < A < 1.5, the preferred values being 9 = 1.4 
(Q > 0.6) and A = 1.3 (Q > 0.98). Our estimates for the 
fitting parameters are = 1.156 ± 0.015, ly = 1.8 ± 0.2, 
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FIG. 2. Finite-size-scaling plot with the form in Eq.(^) for 
a) O = xsG and b) O = The insets represent the same 
data in a log- log plot, showing power-law decay for large ^/L. 



^ = -0.26 ± 0.04, 0.7 ± 0.2, = 0.5 ± 0.3, 

h — 3.3 ± 0.3 and d = 0.9 ± 0.1, where the errors take 
into account the uncertainties on B and A. We then ob- 
tain 7 = v(2 — rj) = 4.1 ± 0.5. As shown in Fig. ^ and 
^, corrections to scaling are important for ^ < 10 [p^ . 
Since the fits do not include the analytic corrections to 
scaling, A and 6 should be regarded as "effective" ex- 
ponents. For comparison, we quote some estimates from 
other MC works: = 1.175 ± 0.025 [|, 1.11 ± 0.04 
[i, 1.13 ±0.06 il, 1.19 ±0.01 |17); v =_1.3 ± 0.1 [|, 



1.20 ±0.04 [|, 1.7 ±0.3 [|, 2.00 ±0.15 §, 1.33 ± 0.05 
r/ = -0.22 ± 0.05 §, -0.35 ± 0.05 §, -0.30 ± 0.06 
-0.37 ±0.04 g, -0.22 ±0.02 Q (notice that in Ref. 
a gaussian distribution of the bonds was considered), 
(ii) We fit our data to 



(8) 



testing the fit stability as above. The fits are excel- 
lent with ^rn > 1-3 but, due to strong correlations be- 
tween a and Tc, the errors on the fit parameters are 
large. For = 1.9 the best fit gives a = 0.20 ± 0.05, 
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Tc = 1.13 ± 0.02, = (1.0 ± 0.2) X 10-3, = 7 ±2 
{Q = 0.77). Notice, however, that any power-law can be 
approximated by an exponential with sufficiently small 
a. For = 3.8 the best fit (shown in Fig. |^) gives 
CT 0.5 ± 0.3, Tc = 1.08 ± 0.04, = (1.1 ± 0.8) x 10-\ 
= 2.4 ± 1.5 (g = 0.69). The deviations of the data 
from this fit for ^ < 3 are consistent with corrections to 
scaling of « 10%. In general, in the presence of an expo- 
nential singularity we expect multiplicative logarithmic 
corrections to Eq.(|^). Our data fit well to 

xsG(o = ^^;e'-"' (logcr (9) 

for ^rn > 2, giving 5/ = 1.30 ± 0.03, ^ = -0.36 ± 0.03, 
r = -0.36 ± 0.06 (Q > 0.9) (see also Fig. |). 
(iii) When we fit our data to 



C(T) = exp(5e/r° 



(10) 



we find that a increases continuously with , from cr « 3 
to cr « 9 |2^ . Even assuming that a stabilizes for higher 
^, we believe that a value ct > 9 is implausibly large. 
In fact, Eq. (|l^) implies a renormalization group (RG) 
transformation dT/dl oc T'^+^ (e' being the RG scale 
factor), while for T — > (at the lower critical dimension) 
we expect dT/dl = a2T'^ + a^T^ -h ... (02 = in the 
phenomenological RG theory of Ref . |^ ) . 

To conclude, we have shown that FSS is verified in 
the 3D Ising spin glass and that the correlation length 
diverges at a finite temperature. Whether this is a con- 
ventional continuous phase transition (in which case the 
lower critical dimension is probably close to three) or a 
transition to a line of critical points, is still not known. 




FIG. 3. Critical behavior of the infinite volume data. The 
solid line is the best fit to Eq.(^ for ^ > 1.8, the dotted line 
is the leading term from the same fit, the dot-dashed line is 
the best fit to Eq.(§l) for ^ > 2.2. 
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FIG. 4. Critical behavior of the infinite volume data. The 
solid line is the best fit to Eq.(^ for ^ > 1.9, the dotted line 
is the leading term from the same fit, the dashed line is the 
best fit to Eq.(^) for ^ > 3.8, the dot-dashed line is the best 
fit to Eq.(|l|) for ^ > 14. 
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